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Abstract 

We investigate the nonlinear power spectra of density perturbations and acoustic oscillations in 
growing neutrino quintessence. In this scenario, the neutrino mass has a strong dependence on the 
quintessence field. The induced coupling stops the evolution of the field when the neutrinos become 
nonrelativistic, and triggers the transition to the accelerating phase of the cosmological expansion. 
At redshifts around five, the neutrino fluctuations are still linear and acoustic oscillations are present 
in the neutrino power spectrum, induced by the acoustic oscillations in the baryonic and dark-matter 
sectors. The neutrino perturbations become nonlinear at redshifts around three. The mode cou- 
pling generated by the nonlinearities erases the oscillations in the neutrino spectrum at some redshift 
above two. There is a potential danger that at later times the influence of the gravitational po- 
tentials induced by the neutrino inhomogeneities could erase the oscillations from the baryonic and 
dark-matter spectra, making the scenario incompatible with observations. For the scenario to be 
viable, the neutrino-induced gravitational potentials in the range of baryonic acoustic oscillations 
should not grow to average values much larger than 10 -4 . The magnitude of the expected potentials 
is still not known reliably, as the process of structure formation is poorly understood in growing 
neutrino quintessence. 



PACS numbers: 95.36. +x, 13.15.+g, 98.80.Es 



1 Introduction 



A popular extension of the quintessence scenario [I] assumes the presence of a coupling between the 
dark-energy and dark- matter sectors [2]. This assumption provides an extended framework in which 
one may hope to address the coincidence problem, i.e. the reason behind the comparable present 
contributions from dark matter and dark energy to the total energy density. In a variation of this 
scenario, characterized as growing neutrino quintessence, the interaction with dark energy is shifted 
from the dark matter to the cosmological neutrino sector [3J. The neutrino-dark energy coupling f3 
can be large, and generate a force substantially stronger than the standard gravitational interaction 
[4]. As a result, even if the neutrinos contribute only a small fraction to the total energy density, they 
may have a significant effect on the cosmological evolution [3] . Their effect becomes important when 
the neutrino mass stops being negligible and the neutrinos become nonrelativistic. This happens at 
a redshift z nr ~ 5 — 10, with the exact value depending on the particular model. At lower redshifts, 
the presence of the neutrinos forces the quintessence field to stop evolving, so that its potential acts 
as an effective cosmological constant, whose value is related to the present neutrino mass. 

The presence of an additional force in the neutrino sector, which is stronger than gravity by 
a factor ~ /3 2 , has profound implications for the evolution of cosmological neutrino perturbations. 
As soon as neutrinos become nonrelativistic (z < z nr ), the perturbations in their energy density 
start to grow, with a characteristic timescale that is shorter by a factor ~ f3~ 2 relative to the one 
characterizing standard gravitational collapse. Static solutions of Einstein's equations are known, 
which describe neutrino lumps held together by the force mediated by the quintessence field. Such 
configurations may be the endpoint of the collapse process 0. The analysis of the growth of neutrino 
perturbations beyond linear level is complicated and the final stages are not well understood [BJ. 
The gravitational potential induced by large-scale structures in the neutrino sector, which can have 
a size of 100 Mpc or more, can affect the propagation of photons and influence the cosmic microwave 
background (CMB) through the integrated Sachs-Wolfe (ISW) effect. An analysis based on linear 
perturbation theory leads to the conclusion that models with a strong neutrino-dark energy coupling 
are excluded [7j . However, an extrapolation of linear growth is strongly misleading, as it would imply 
completely unrealistic neutrino-induced cosmological gravitational potentials 8 . In particular, the 
effects of backreaction and virialization are not accounted for in a linear treatment. If these are 
sufficiently strong, the neutrino structures may not be as dense as the extrapolation of the linear 
analysis would indicate. The resulting gravitational potentials may be sufficiently small so as not to 
affect the CMB significantly [8]. As a general rule of thumb, potentials larger than ~ 10 at length 
scales of 100 Mpc or more give too strong an effect on the CMB. 

In this work we study the power spectra of dark-matter and neutrino perturbations in the scenario 
of growing neutrino quintessence. The growth of perturbations can be used in order to constrain the 
scenario through comparison with the observed large-scale structure. The most promiment feature 
of the baryonic and dark-matter spectrum is a series of peaks and valleys, characterized as baryon 
acoustic oscillations (BAO). They originate in the period of recombination, and correspond to sound 
waves in the relativistic plasma of that epoch. The chacteristic length scale of BAO is around 100 
Mpc. Even in the standard ACDM scenario, the exact form of the dark-matter power spectrum 
at such scales is not easy to compute precisely, because of the failure of linear perturbation theory 
to describe reliably the growth of the corresponding fluctuations under gravitational collapse. At 
length scales below about 10 Mpc, the evolution is highly non-linear, so that only numerical N-body 
simulations can capture the dynamics of the formation of galaxies and clusters of galaxies. However, 
fluctuations with length scales of around 100 Mpc fall within the mildly non-linear regime, for which 
analytical methods have been developed. We focus on scales in the range 50-200 Mpc, within which 
BAO are visible. In growing neutrino quintessence the neutrino power spectrum displays a much 
faster growth and overtakes the dark-matter spectrum at redshifts below z ~ 4. The nonlinear 
corrections become very large, so that analytical methods become unreliable even in the BAO range 
at redshifts near z = 0. Our aim is to explore the range of validity of the analytical methods and 
investigate the form of the spectra within this range. 

The various analytical methods [9]-[l6] that have been developed in order to go beyond linear 
perturbation theory essentially amount to resummations of subsets of perturbative diagrams of arbi- 
trarily high order, in a way analogous to the renormalization group (RG). We follow the approach of 
[13] . named time-RG or TRG, which uses time as the flow parameter that describes the evolution of 
physical quantities, such as the spectra. The method has been applied to ACDM and quintessence 
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cosmologies [13], allowing for a possible coupling of dark energy to dark matter |17| . or a variable 
equation of state |18| . It has also been used for the study of models with neutrinos of constant mass 
[19] . A comparative analysis of several analytical methods, using N-body simulations as a reference, 
has been carried out in ref. [5D] for ACDM cosmology. The study demonstrates that TRG remains 
accurate at the 1-2% level over the whole BAO range at all redshifts. 

In the following section we summarize the formalism we use. In section 3 we present the results 
of the numerical integration of the evolution equations for the spectra. Finally, in section 4 we give 
our conclusions. 

2 Neutrinos coupled to dark energy 
2.1 Evolution equations for the perturbations 

We assume that the energy density of the Universe receives significant contributions from three 
components: a) nonrelativistic matter, in which we group standard baryonic matter (BM) and cold 
dark matter (CDM); b) massive neutrinos, whose equation of state varies during the cosmological 
evolution; c) a slowly varying, classical scalar field 0, named cosmon pQ, whose contribution to the 
energy density is characterized as dark energy (DE). We also allow for a direct coupling between 
the neutrinos and the cosmon field. In particular, we assume that the mass m„ of the neutrinos 
has a dependence on <j>, and define the coupling as P{4>) = —dlnm v (<f>)/d(f). In the present paper we 
consider a model with constant f3. The equation of motion of the cosmon field takes the form 



At early times, when the neutrinos are relativistic and their energy-momentum tensor is traceless, 
they are effectively decoupled from the cosmon field. Only when the neutrinos become nonrelativistic 
the cosmon-neutrino coupling f3 becomes effective, leading to energy exchange between the neutrino 
and dark-energy sectors. We normalize all dimensionful quantities, such as the cosmon field, with 
respect to the reduced Planck mass M = (87rG) _1/ ^ 2 . This is equivalent to setting M = 1. 

We are interested in the more recent stages of the cosmological evolution (redshifts z 10). In 
the scenario we consider, the neutrinos become nonrelativistic at such redshifts. In the remaining of 
the section, in which we develop the formalism for the treatment of the nonlinear corrections to the 
power spectra, we assume that the neutrino pressure vanishes. We approximate the metric as 



We assume that the Newtonian potential $ is weak: $ < 1. Also we decompose the cosmon field 
4> as 4>(t,x) — 4>(t) + 84>(t,x), with 8<j> <C 1. In general, (j> = 0(1) in units of M. Finally, we 
decompose the density as p{r,x) — p(r) + Sp(r,x), while we neglect the pressure. A self-consistent 
expansion scheme can be obtained if we assume the hierarchy of scales: $, 5cj> <C \Sv\ <C Sp/p < 1. 
Such a hierarchy is predicted by the linear analysis for subhorizon perturbations with momenta 
k 2> H — a/ a. For subhorizon perturbations, it is consistent to make the additional assumption 
that the spatial derivatives of $, 5<j) dominate over the time derivatives. The predictions of the linear 
analysis allow us to make a more quantitative statement. We assume that a spatial derivative acting 
on <3>, 5<j} or Sv increases the position of that quantity in the hierarchy by one level. In this sense V$ 
is comparable to H8v, while V 2 $ is comparable to H 2 5p/p. 

With the above assumptions, one can derive the equations that describe the evolution of the Uni- 
verse. The details of the calculation are given in [17] • The evolution of the homogeneous background 
is described by 
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where we have defined ptot = T^iPi + <P/(2ci 2 ) + U(<f>). The values i = 1, 2 correspond to neutrinos 
and CDM+BM, respectively. For the neutrinos we consider a constant coupling f3\ — /3, while we 
assume that the CDM+ BM sector does not couple to the cosmon field: /3 2 = 0. 

We describe the perturbations in terms of the cosmon perturbation 5(f), the Newtonian potential 
<3>, the density perturbations 8 pi and the velocity fields Vi. We have two Poisson equations 



v 2 50 = - a 2 fcs Pl = -3 p t n 2 n t s 1 

i i 



(2.4a) 
(2.4b) 



with Qi(r) = pid 2 /(3H 2 ), and the continuity and Euler equations 

8p d + SHSpi + V[(pi + 5pi)8vi] = — Pi<j>5pi 
8% + (U- Pi$)Si!i + (Si>i ■ V)5vi = - V$ + PiS76<t>. 



(2.5a) 
(2.5b) 



The inspection of eqs. (|2.4ap - (12.5b[) reveals a potential shortcoming of the assumed hierarchy 
for large f3. From the Euler equation one concludes that the scalar quantity comparable to the 
gravitational potential is actually /36<j>. From eqs. (|2.4ap . (|2.4b|) we can see that (58(f) is typically 
larger than $ by a factor ~ /3 2 . This limits the range of validity of a hierarchy in which /38(f) and $ 



are treated on equal footing. Furthermore, eqs. (|2.5a[) . (|2.5b|l implicitly neglect effects arising from 
a difference between the local neutrino mass, as given by the local value of the cosmon field, and the 
average cosmological neutrino mass, as given by the background value <j>. Including higher orders in 
wold increase considerably the algebraic complexity of the evolution equations. 



2.2 Evolution equations for the power spectra 

The evolution equations are expressed in their most useful form in terms of the density contrasts 
Si = Spi/pi < 1 and the divergence of the velocity field 9i(k,r) = V ■ 8vi(k, r). For the Fourier 
transformed quantities, we obtain from eq. (|2.5a|l 
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Eqs. (f23b|) . Igga]! , (f2~4b|) give 
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where r\ — ln[a(r)/oi], with a,i the scale factor at some convenient time, at which we define the initial 
conditions. (We normalize the scale factor so that ao = 1 today.) The index v refers to neutrinos 
and m to the CDM+BM sector. This allows us to bring eqs. Ipl!)). (f278|) in the form p [T2l[T3] 

<9,,V?a(k, 77) + Q ab ip b (k, 77) = e' 7 7 tl i )c (k, -ki,-k 2 )^6(ki,?7)(/p c (k;2,?7). (2.11) 

The indices a,b,c take values 1,...,4. The values 1, 2 characterize neutrino density and velocity 
perturbations, while 3, 4 refer to CDM+BM quantities. Repeated momenta are integrated over, 
while repeated indices are summed over. The functions 7, that determine effective vertices, are 
analogous to those employed in |12U13| . The non-zero components are 



7121 (k, ki,k 2 



a(ki,k 2 



-fo(k + ki +k 2 ) = 7 ii 2 (k,k 2 ,ki) 



The il-matrix is 



7 222 (k,ki,k 2 ) = /3(ki,k 2 ) 5 D (k + ki + k 2 ) 

7343(k, k 3 , k 4 ) = " ^' ^ 5d (k + k 3 + k 4 ) = 7334 (k, k 4 , k 3 ) 



7444(k, k 3 ,k 4 ) = /3(k 3 ,k 4 ) 5 D (k + k 3 + k 4 ) 
/ 1 



(2.12) 



-1 



-0,(2/3 2 + l) 2-/30' + 



H 



V 











- f2. n 



7 f^T7 



2 + 






-1 
H' 



(2.13) 



H 



) 



where a prime denotes a derivative with respect to rj. 

The next step is to derive evolution equations for the power spectra. The spectra and bispectra 
are defined as 



(f>a(K r))cp b {q, 77)) =5 D (k + q)P a t(k, rj) 
{<p a {k, T))tp b (q, T))ip c (p, 7])) =5 D (k + q + p)B abc (k 1 q, p, rj). 



(2.14) 



The essential approximation that we have to make in order to obtain a closed system of equations is 
a truncation of the four-point function which appears in the evolution equation for the trispectrum. 
In this way we obtain [131 117] 

d v P ab (k, 77) = -Q ac P cb (k, 77) - D. bc P ac (k, 77) 

+e v J d 3 q[y acd (k, -q, q - k)B bcd (k, -q, q - k) 

+7i>cd(k, -q, q - k)B acd (k, -q, q — k)] , (2-15) 
d v B abc (k., -q, q - k) = -Q ad B db c(k, -q, q - k) - f2 6d _B adc (k, -q, q - k) - Q, cd B ab d{\i, -q, q - k) 

+2e 7) / d 3 g[7 ade (k, -q, q - k)P db (q, 77)P ec (k - q, 77) 



(2-16) 



+7Me(-q,q- k,k)P dc (k- q,77)P ea (k, 77) 
+7cd e (q - k, k, -q)P da (k, 77)P ei ,(q, 77)] . 



3 Numerical analysis 
3.1 Approximations 

The presence of two massive species (neutrinos and CDM+BM) complicates the structure of the 
equations compared to the case where they are treated as a single fluid, discussed in [T3]. The full 
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Figure 1: The fractional energy density in neutrinos (solid), CDM+BM (dashed) and dark energy (dot- 
ted). The acceleration parameter q — ad /a 2 (dot-dashed) is also depicted. 



system of eqs. (|2.15p . (|2. 16p contains 74 equations, namely, 10 for the power spectra and 64 for the 
bispectra, compared to the 11 equations of the single-matter case [13]. An accurate calculation also 
requires the discretization of the fe-space with at least 500 points. Taking into account that the 
bispectra depend on three external momenta, it is apparent that the necessary computing power is 
significant. In order to make the numerical integration of the evolution equations feasible, we have 
to make additional approximations. 

As will be apparent in the following, the method we are employing remains valid only within 
the mildly nonlinear regime. Strong nonlinearities and backreaction effects, as well as the process of 
virialization, are not accounted for in our treatment. In the scenario we are considering, the strongest 
nonlinearities appear in the neutrino sector, while the CDM+BM spectrum remains within the linear 
regime. It is, therefore, a good approximation to neglect the nonlinear terms in the evolution equations 
for the CDM+BM spectrum. In practice, this means that we may set the vertices 7334, 7343, 7444, 
defined in eq. (|2.12p . equal to zero. 

We derive an approximate solution of eqs. (|2.15[) . (|2.16|l in the following way: 

• At a first stage, we integrate the full system of equations setting all the vertices y a tc equal to 
zero. This reproduces the linear spectra. 

• Subsequently, we integrate the 11 equations for P a t, B a t c with a,b = 1,2. Notice that these 
involve only the vertices 7112, 7121, 7222- All spectra and bispectra appearing in these equations 
with an index 3 or 4, because of the nonzero entries ^23 and ^41 entries of eq. (|2.13[1 . are 
approximated by their linear solutions derived at the previous stage. 

The above procedure can reproduce the modifications in the neutrino spectrum induced by mode 
coupling at the nonlinear level. Our main result will be the effect of this process on the BAO in 
the spectrum. Our analysis is performed at redshifts above 2 ~ 2, at which the only significant 
nonlinearities appear in the neutrino sector. 

3.2 Results 

We consider a model in which the cosmon field has a potential V(<j)) = M 4 exp(— with a = 10. As 
we mentioned earlier, we normalize dimensionful quantities with respect to the Planck mass M = 1. 
The coupling between the neutrinos and the cosmon field is taken to be constant. This is equivalent 
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Figure 2: The neutrino density power spectrum Pn(k,rf) (solid lines) and the CDM+BM density spec- 
trum P 33 (k,r)) (dotted lines) at redshifts z = 4.70, 4.08, 3.04, 2.77, 2.69, 2.60 (starting from below). 

to assuming that the neutrino mass is m„((/>) = m„ exp(— j3(j>). In our calculation we consider three 
degenerate neutrinos and use /3 = —52. The constant rh v and the present value of the background field 
4> are chosen such that the present neutrino mass is m„o = 2.1 eV. The evolution of the cosmological 
background for z < z nr = 8 is shown in fig. [T] We have approximated the neutrinos as nonrelativistic 
during this period. We depict the fractional energy density in neutrinos (solid), CDM+BM (dashed) 
and dark energy (dotted). We also plot the acceleration parameter q = ad/ a 2 (dot-dashed). The 
cosmological expansion becomes accelerating at a redshift z ~ 1. The present fractional contributions 
of neutrinos, CDM+BM and the cosmon field to the total energy density are f2„ = 0.13, f2 m = 0.23, 

= 0.64, respectively. The spatial curvature is assumed to vanish, while the current expanstion 
rate is H = 72 km s _1 Mpc -1 . 

We start the integration of the evolution equations for the spectra at a redshift z nr = 8. The 
neutrinos are already nonrelativistic at this time in our model. The initial conditions at z nr are 
obtained through the implementation of the background and linear-perturbation equations in the 
code CAMB [21], generalized for the case of mass- varying neutrinos. We assume that the primordial 
spectrum is scale invariant with spectral index close to 1. In the background evolution, we take into 
account the transition of neutrinos from being relativistic at high redshifts to nonrelativistic near z nr . 
We point out that the neutrinos are essentially decoupled from the cosmon field during the time that 
they are relativistic (see eq. (12. l[) b We perform a careful analysis of the evolution of the neutrino 
perturbations by solving the Boltzmann equation. In this way we obtain the neutrino and matter 
spectra at z nr . The mixed spectra are obtained as the geometrical averages of the pure ones, i.e. 



etc, consistently with the expectation in the linear regime. 

The evolution of the spectra at redshifts below z nr is linear to a good approximation down to 
z ~ 3. Around this time the neutrino spectrum has grown sufficiently for the nonlinear corrections 
to affect the evolution. This is demonstrated in fig. [2] where we depict the neutrino density power 
spectrum Pn(k, rj) (solid lines) and the CDM+BM density spectrum P33(fc, 77) (dotted lines) at various 
redshifts. The solid lines, starting from below, correspond to z = 4.70, 4.08, 3.04, 2.77, 2.69, 2.60. 
The dotted lines correspond to the same redshifts, but they almost coincide. We point out that the 
actual density spectra are obtained from those depicted in fig. [2] by multiplication with a factor 
exp(27?) = (a(z) / a(zi)) 2 = 81(1 + z)~ 2 , where we have taken Zi = z nr — 8 in our calculation. This 
factor arises because of the definitions (|2,10[) . (|2.14[) . while 77 = \n[a(r)/ai], with a; the scale factor 
at some convenient time, at which we define the initial conditions. In our calculation a,i — 1/9. The 
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spectra are given in units of (Mpc//i) 3 , consistently with their definition (|2.14p . They are also smaller 
by a factor (2-7r) 3 compared to the definition in CAMB. 

At high redshifts the neutrino density spectrum is very suppressed because of free-streaming, as 
the neutrinos are relativistic over the entire depicted momentum range. It is apparent from fig. [5] 
that, as soon as the neutrinos become nonrelativistic, their spectrum grows very fast because of the 
additional attractive force generated by the neutrino-dark energy coupling. This force is 2/3 2 ~ 5400 
times stronger than gravity. The neutrino density spectrum at z = 4.70, 4.08, 3.04 is within the 
linear regime, while the nonlinearities become apparent at z —2.77, 2.69, 2.60: starting from large 
k, the neutrino spectrum grows faster than the linear analysis would predict. The matter spectrum 
retains its shape, as its evolution is described by linear theory to a very good approximation for the 
whole range of redshifts depicted in fig. [2] 

The most interesting feature of the spectra in fig. [5] is their oscillatory behavior. The BAO in 
the matter sector originate in the period of recombination, and correspond to sound waves in the 
relativistic plasma of that epoch. The corresponding oscillations in the neutrino spectrum develop as 
soon as the neutrinos become nonrelativistic and start falling in the gravitational potentials generated 
by the CDM+BM inhomogeneities. The corresponding process takes place within the linear regime 
and is fast, as the CDM+BM power spectrum exceeds that of neutrinos by several orders of magnitude. 
The amplitude of oscillations in the neutrino spectrum starts diminishing as soon as the nonlinear 
corrections become important. Again the transition is fast, and the oscillations are hardly visible at 
z = 2.60. This phenomenon has its origin in the mode coupling induced by the nonlinearities. It 
is observed within the ACDM model as well: near z = the higher peaks in the matter spectrum 
disappear and only the first two or three survive when the nonlinear corrections are taken into account. 
In the model at hand, the nonlinearities in the neutrino sector are very strong, because of the rapid 
growth of the spectrum. The oscillations disappear at a redshift above z ~ 2. 

A crucial question is whether the disappearance of oscillations from the neutrino sector induces 
their elimination from the CDM+BM sector as well. We address this issue in the following section. 
An important quantity in the context of this discussion is the gravitational potential induced by the 
inhomogeneities. We can have an estimate of its magnitude at various scales by starting from the 
Poisson equation (|2.4b[) for the gravitational potential and defining an associated power spectrum as 

A|(fe,r?) = 47vk s P*(k,r,) = 47rfc 3 e 2 " (~^) (njfti + 2Q. v Q, m P 12 + H^P 22 ) . (3.2) 

In fig. [3] we plot A* as a function of k for redshifts z = 4.70, 4.08, 3.04, 2.77, 2.69, 2.60. At 
2 = 4.70 the potential is mainly generated by the CDM+BM inhomogeneities, while for z = 2.60 
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the main contribution at large k comes from the neutrino inhomogeneities. This is apparent in the 
change of shape of the function A$(fc). Already at this stage the potential takes values in the region 
[1CT 6 - lCT 5 ] in the BAO range. At z 2.6 the neutrinos dominate the gravitational potential, so 
that its magnitude depends on the growth of inhomogeneities in this sector. We point out that k is 
the comoving momentum, so that the physical length scale is related to k/a(rf). 

The form of the spectra at redshifts below z ~ 2.60 cannot be obtained within the scheme we are 
employing, as the validity of TRG lies within the mildly nonlinear regime. At its present development, 
the scheme cannot account for the process of virialization and formation of bound structures. As a 
result, inhomogeneities tend to grow without limit, and the power spectrum diverges. The integration 
of eqs. (I2.15p . (|2.16p below z = 2.60 results in an increase of the spectrum much faster than the one 
predicted by the linear treatment. 

4 Conclusions 

The numerical analysis of the previous section has led to some concrete conclusions: 

• The neutrino power spectrum is subdominant to the CDM+BM spectrum at high redshifts 
because of free-streaming during the period that neutrinos are relativistic. The neutrino mass 
grows during the cosmological evolution of the cosmon field. As soon as the neutrinos become 
nonrelativistic (z ~ 8 in our model), their spectrum grows rapidly and overtakes the CDM+BM 
spectrum. 

• During the period that the neutrinos are nonrelativistic and the evolution linear (8 <; z <; 3), 
the neutrino spectrum develops oscillations induced by the BAO in the CDM+BM sector. The 
oscillations in the neutrino sector are erased by the mode coupling generated by the nonlinearities 
when these become significant (3 <; z <; 2.6). 

• At high redshifts (z <; 3) the gravitational potentials in the BAO range are dominated by 
CDM+BM inhomogeneities. At lower redshifts, the neutrino inhomogeneities start dominating 
and they determine the potentials. At the time when the perturbations in the neutrino sector 
become highly nonlinear [z ~ 2.6 in our model), the potentials take values in the region [10 -6 — 
10" 5 ] in the BAO range. 

The issue that cannot be resolved by our analysis is whether the BAO persist in the CDM+BM 
sector at redshifts below the one at which the neutrino sector becomes highly nonlinear. An ex- 
trapolation of our results below z = 2.6 within a linear treatment would indicate that the potentials 
induced by the neutrino inhomogeneities grow rapidly extremely large, so that their influence would 
erase the BAO from the CDM+BM sector already at redshifts above z = 2. Such a conclusion would 
make growing neutrino quintessence incompatible with observations. However, the linear increase is 
strongly misleading, as the neutrino-induced gravitational potential would exceed the extreme value 
resulting from all neutrinos within the horizon being concentrated at a single point [§]. Unfortunately, 
our present approximation cannot account for a slowing down of the growth as compared to the linear 
approximation. 

Nevertheless, an important lesson can be learned from a linear extrapolation. One finds that, 
at the time when the oscillatory behavior disappears from the CDM+BM density spectrum, the 
gravitational potentials in the BAO range take values in the region [10~ 3 - 10~ 2 ]. A significant 
reduction of the amplitude of oscillations, at the 10-20% level, corresponds to potentials in the region 
[10" 4 -10~ 3 ]. The same critical orders of magnitude of the gravitational potential persist even if the 
neutrino spectrum is modelled artificially to grow at a much slower rate than the one predicted by 
an extrapolation of our results: the disappearence of BAO at some redshift between z ~ 2 and z ~ 
is linked to the presence of gravitational potentials ~ 10~ 3 or larger. It has been established already 
that the appearance of such potentials would make the predictions of growing neutrino quintessence 
for the CMB incompatible with its measured spectrum, because of a very strong ISW effect [8]. 

We expect our qualitative conclusions to remain valid for all variations of the scenario of growing 
neutrino quintessence (different forms of the potential U(<j>) and the coupling P(<f>))- The strong 
influence of the neutrino sector on the background evolution requires a large coupling to the cosmon 
field. In turn, this generates a strong long-range force between neutrinos. This causes the neutrino 
perturbations to become nonlinear early in the cosmological evolution, and the oscillations to disap- 
pear from their density power spectrum. The influence on the CDM+BM sector and the observable 
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BAO depends on the final size of nonlinear inhomogeneities in the neutrino sector and the respective 
gravitational potentials. Similar considerations would also apply to other models of mass-varying 
neutrinos [4], if the neutrinos influence the cosmological evolution considerably. 

Our analysis of density spectra in growing neutrino quintessence sets a rough quantitative bound 
on the typical gravitational potentials generated by the neutrino inhomogeneities: they should not 
become much larger than 10 -4 in the BAO region. Otherwise, the amplitude of BAO will be reduced 
or completely eliminated. As has been discussed in [5J|5], the formation of structures in growing 
neutrino quintessence is not well understood. It is possible that highly nonlinear structures or virial- 
ized objects start forming early, and their backreaction prevents the neutrino power spectrum from 
growing very large. It is conceivable that the gravitational potentials at large scales remain suffi- 
ciently small for the CMB spectrum and the BAO to remain largely unaffected. However, a better 
understanding of structure formation in the scenario of growing neutrino quintessence is required in 
order to resolve this issue. 
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